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Abstract 

We present the result of a search for charged Higgs decays of the top quark, produced in pp cohisions at -y/s = 
1.8 TeV. When the charged Higgs is heavy and decays to a tau lepton, which subsequently decays hadronically, the 
resulting events have a unique signature: large missing transverse energy and the low-charged-multiplicity tau. Data 
collected in the period 1992-1993 at the Collider Detector at Fermilab, corresponding to 18.7±0.7 pb~^, exclude 
new regions of combined top quark and charged Higgs mass, in extensions to the standard model with two Higgs 
doublets. 



I. INTRODUCTION 

We have conducted a search for decays of the top 
quark to a charged Higgs boson, using the Higgs de- 
cays to hadronically decaying tau leptons. The re- 
sults presented here come from data collected during 
the years 1992-1993 at the Collider Detector at Fer- 
milab, corresponding to an integrated luminosity of 
18.7±0.7 pb~^. A charged Higgs arises in extensions 
to the Standard Model with two Higgs doublets 
If the charged Higgs exists in such a model and is 
lighter than the top quark, then two competing chan- 
nels are possible: t^H'^h and t^W~^h. The charged 
Higgs can decay either to tu or to cs. The branching 
ratios of these processes depend on the top quark 
and charged Higgs masses, and on tan /?, the ratio 
of the vacuum expectation values of the two Higgs 
doublets in the modelJil We consider here only the 
kinematically allowed cases where m-top > mn + rrih 
and mtop > mw + n^b- In these cases three decay 
modes of the quark pairs are possible: tt^H^H^bb, 
ti-^H^W'^bb, and tt^W+W-bb. 

This analysis of hadronic decays of the r-lepton 
uses a method similar to that used previously [§] but 

^We consider here only models in which one Higgs doublet 
couples to the up-type quarks, and the other doublet couples 
to the down-type quarks and the leptons. 



with four times more integrated luminosity and an 
event selection designed for larger top and charged 
Higgs masses. Our most recent limit uses the 
leptonic decays of the tau using the data set col- 
lected in the same period as in this paper, i.e. 1992- 
1993. Charged Higgs masses from 45 GeV/c^ to 110 
GeV/c^, and top quark masses from 90 GeV/c^ to 
110 GeV/c^ were excluded at a 95% CL, as shown in 
the lower hatched part of the plot in Fig. ||. Exper- 
iments at LEP exclude a charged Higgs with mass 
less than 45 GeV/c^ §. 

In the present analysis a more stringent limit 
results from the 64% hadronic branching ratio of 
the tau, compared with the 36% leptonic branch- 
ing ratio. However, the larger expected back- 
ground, mainly hadronic processes, must be well 
modelled. The analysis presented here addresses 
top masses in the range extending from the lim- 
its of previous searches ^, about 100 GeV/c^, 
up to the mass range which has been measured, 
176±8(stat)±10(syst) GeV/c^ g. The analysis ex- 
cludes by direct search a top or top- like object de- 
caying via a charged Higgs in this region. 

Top quark pair events with one or two charged 
Higgs decays should contain energetic jets which 
come from b quarks and the decays of the taus. 
Each top quark leads to the production of two en- 
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ergetic neutrinos, leading in turn to a large missing 
transverse energy, denoted The first neutrino 

is emitted in the charged Higgs or W boson decay, 
and the second results from the tau decay. The sig- 
nature of an hadronically decaying tau is a narrow 
jet with either one or three associated charged parti- 
cles. Thus, top quark pair events with charged Higgs 
decays to tau leptons can be found by looking for 
an excess of events with narrow energetic jets with 
one or three charged particles coming from the tau, 
along with the presence of other energetic hadronic 
jets and neutrinos. This signature differs from that 
of standard model top events due to the higher prob- 
ability of decays to tau leptons, and larger missing 
Ex- 

II. CDF DETECTOR AND TRIGGER 

The CDF detector is described in detail else- 
where The most important components of the 
CDF detector for this analysis are the tracking cham- 
bers and calorimeters. The relevant tracking cham- 
bers are the vertex time projection chamber (VTX) 
and the central tracking chamber (CTC), which is 
a large cylindrical drift chamber surrounding the 
VTX. Both are located inside a superconducting 
solenoid magnet generating a 1.4 T field. The VTX 
provides z-vertex reconstruction and r — z tracking 
over the pseudorapidity range \r]\ < 3.25 [0], where 
the z axis is the proton direction along the beam 
line and r refers to the radial coordinate transverse 
to the beam line. The momenta of charged parti- 
cles are measured in the CTC. The solenoid and the 
tracking volume of CDF lie inside electromagnetic 
and hadronic calorimeters which cover 27r in azimuth 
and up to \ri\ = 4.2. The calorimeters are segmented 
in azimuth and pseudorapidity to form a tower ge- 
ometry which points back to the nominal interaction 
point z = 0. 

The "trigger" decision as to whether or not the 
data from a particular interaction should be recorded 
depends on the particular pattern of energy de- 
posited in the calorimeters, the presence of charged 
tracks in the CTC, and the presence of penetrating 

^For a calorimeter energy deposit, assuming the particles 
came from some point along the beam axis, a direction in 
space is defined. The transverse energy Et is the component 
of the energy vector in the plane perpendicular to the beam 
axis. The is defined as the magnitude of the vector sum of 
the transverse energy Et of each calorimeter energy deposit. 



charged particles in the muon chambers which sur- 
round the calorimeter. This analysis relies in par- 
ticular on a trigger which uses analog sums of the 
calorimeter energy deposits to determine the missing 
transverse energy. Since the charged Higgs events 
sought in this analysis generally have large miss- 
ing transverse energy, this analysis uses only those 
events which satisfy a trigger requirement of at least 
35 GeV of missing transverse energy. 

HI. EVENT SELECTION 

The criteria to reject background and to select the 
charged Higgs signal were determined using a Monte 
Carlo simulation based on top quark and charged 
Higgs masses just beyond those excluded in pre- 
vious analyses, namely mtop = 120 GeV/c^, and 
"^Higgs = 100 GeV/c^. A version of the Monte Carlo 
program IS A JET Q, modified to correctly model 
the polarization of the taus, generated events which 
were then passed through the CDF detector simula- 
tion. 

The selection criteria aim to select events with 
large missing transverse energy due to the neutrinos, 
the presence of a hadronically decaying tau lepton, 
and at least one other jet due either to another tau 
or to one of the jets from the top quark. Each event 
must have 

• > 40 GeV, 

• 5(^t) > 4 GeV^/2^ 

• a tau lepton identified as discussed below, with 

- -Et > 30 GeV, 

- < 1, 

• a jet as defined below, with 

- Et> 20 GeV, 

- 1^1 <2, 

• A4>r-jet < 140°, and 

• scalar T.\Et\ > WO GeV, 

where we use the definition S{^t) = ^t/ \/'^\Et\ 
for the "significance" of the missing Ej^. 

The number of events which satisfy each criterion 
are listed in Table I. The relative efficiencies between 
consecutive cuts for the ISAJET Monte Carlo sim- 
ulation with mtop = 120 GeV/c^, and mniggs = 100 
GeV/c^ are also shown in Table I. 
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Figure 1: Distribution of for the data sam- 
ple (solid line) and for the Monte Carlo simulation 
(dashed and dotted line). Events to the left of the 
arrow are removed by the cuts on the data sample. 
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Figure 2: Distribution of S{^t) for the data sam- 
ple (solid line) and for the Monte Carlo simulations 
(dashed and dotted line). Events to the left of the 
arrow are removed by the cut. 



The scalar nature of the charged Higgs implies 
that the two neutrinos produced in the decay chain 
tend to travel in the same direction, resulting in a 
large ^t- Furthermore, the charged Higgs decays 
mainly to a tau for large values of the parameter 
tan p. For smaller tan (3 values, the probability for 
the top quark to decay to a W boson increases, and 
the charged Higgs decays more often to a quark- 
antiquark pair. In this case the average conse- 
quently becomes smaller. Thus, the trigger and se- 
lection requirements on enhance the acceptance 
in the case of large values of tan /3. 

The criteria on the missing transverse energy and 
significance reinforce the trigger requirements and 
select events with energetic neutrinos. The distri- 
butions of ^TjSLiid S{^t) are shown in Fig. ^ and 
Fig. ^, for the data sample and for the Monte Carlo 
simulation of the signal. 

When a top quark decays to a charged Higgs, a 
large fraction of its energy goes into creation of the 
charged Higgs. The smaller remaining energy for 
the b quark produces jets of lower Et- Since the 
charged Higgs carries a large energy, its decay prod- 
ucts receive a strong boost. In particular, the taus 
which come from the charged Higgs have very large 
Et, resulting in a large-transverse- momentum (pt) 
associated charged particle near the jet axis. The 
average pT increases with increasing charged Higgs 



mass. A charged particle is associated with a jet 
if its initial direction points within a cone of radius 
AR = \/ Ar]"^ + A(/)2 = 0.4 of the jet direction, where 
(p is the azimuthal angle in the plane perpendicular 
to the beam axis. 

The event selection requires the presence of at 
least two jets, formed from calorimeter energy de- 
posits in a cone of radius AR = 0.4. The first jet 
must have Et > 30 GeV, lie in the region \7]\ < 1, 
and have an associated charged particle with p-j- > 4 
GeV/c. If more than one jet satisfies these criteria, 
the jet with the largest Et is chosen. 

The requirements on the second jet are less strin- 
gent. The Monte Carlo simulation shows that the 
second jet has a smaller Et, and is less often in the 
central region of the detector, |r/| < 1. The second 
jet must have Et > 20 GeV, \r]\ < 2, and an associ- 
ated charged particle. 

Requiring that the z-intercept of the largest-pr 
charged particle associated with each jet be within 
5 cm of the primary z- vertex of the event rejects jets 
from additional interactions in the event. 

Subsequent criteria to identify hadronically decay- 
ing taus assign one of the two jets as the "tau." Since 
the charged particles in hadronically decaying en- 
ergetic taus must lie in a narrow cone around the 
calorimeter energy deposit, the tau must satisfy an 
isolation criterion in which there must be no asso- 
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Table I: Number of events selected and signal efficiency for each selection criterion. The efficiencies represent 
the successive effect of each criterion, for events in a Monte Carlo simulation with mtop = 120 GeV/c'^ 
and mHiggs = 100 GeV/c^. 



dated charged particle with > 1 GeV/c found 
between cones of 10° and 30° defined around the di- 
rection of the associated track with the largest pT- 
If the ffist jet fails the cut, the algorithm applies the 
cut to the second jet which in addition must then 
pass the stricter jet Et, associated charged particle 
Pt, and rj cuts of the first jet. The jet which passes 
the isolation cut is called the tau candidate. 

A small fraction of electrons and single hadrons 
or low-multiplicity hadronic jets also satisfy the tau 
selection criteria. To reject electrons, the tau can- 
didate must satisfy 1 — {IQEt /'^\pt\ — 1)^^ > Jem-, 
where fEM is the fraction of the total energy de- 
posited in the electromagnetic calorimeter and Ti\pt\ 
is the sum of the magnitudes of the transverse mo- 
menta of charged particles in the 10° cone around the 
jet axis. A fraction of the hadronic background is re- 
jected by a similar cut: 1 — (^E^t/^^Pt — 1)~^ < Jem-, 
where the factor h has been optimized as a func- 
tion of the Et of the tau candidate: b = 0.815 
for 30 GeV< Et < 45 GeV, b = 0.995 for 45 
GeV< Et < 69 GeV, and 6 = 0.860 for Et > 69 
GeV. 

Lastly, the requirement that the scalar sum of 
transverse energy exceeds 100 GeV removes back- 
ground from W -|- jets events in which the W decays 
leptonically. 



IV. TAU SIGNAL IN HADRONIC EVENTS 

In order to demonstrate that the criteria select 
hadronic tau decays in a process known to contain 
tau leptons, one can extract the tau signal from the 
process pp ^ W + jets, W ^ tv hy selecting a sam- 
ple of events with a tau, a jet and missing transverse 
energy. Employing less stringent cuts on the ^t, 
S{^t), Et and Pt of the charged particles associated 
with the jets, and then requiring < 85 GeV 

and a tight cut on the width of the energy deposit 
if the jets, results in the multiplicity distribution 
shown in Fig. |^. The plot shows a clear excess of one 
and three charged particles in the distribution of the 
associated charged particle multiplicity distribution 
in the 10° cone, attributed to W plus jets events. 
The data (points) are compared to a background es- 
timate (cross-hatched histogram), and a HERWIG 
Monte Carlo simulation Q of ^ rzv (open his- 
togram), normalized according to the HERWIG pre- 
dicted cross section. 

V. MONTE CARLO SIMULATION AND 
TRIGGER EFFICIENCY 

The simulation of charged Higgs events uses the 
physics generator ISAJET and the CDF detector 
simulation. In order to compute the acceptance at 
any point in the top quark versus charged Higgs mass 
plane, it is necessary to make a mixture of events 
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Figure 3: The tau signal found in the data sample 
using a similar algorithm with less stringent cuts. 
The data (points) are compared with an estimate of 
jet background (cross-hatched), and with a Monte 
Carlo simulation of W+jets events, with W ^ tv 
(open histogram). 

from the simulation of tt^H^H^bb, tt^H'^W'^bb, 
and tt^W~^W~bb processes for charged Higgs 
masses in the range of 50 GeV/c^-160 GeV/c^ and 
for top quark masses in the range of 100 GeV/c^- 
174 GeV/c^. The tt production cross-section is 
taken from a next-to-leading order theoretical cal- 
culation ijl^ . The simulation of the effect of the 
trigger efficiency comes from a measurement of the 
efficiency as a function of the in events which 
triggered on the presence of jets. 

In the Monte Carlo simulation, more than 95% of 
tau candidates found by the algorithm correspond in 
spatial direction to the actual taus generated in the 
event. 



VI. BACKGROUNDS 

The dominant backgrounds to charged Higgs 
events in the selected event sample are hadronic pro- 
cesses, and processes in which a Z or W is produced, 
possibly accompanied by jets. In almost all of the 
background, an hadronic jet fluctuates to have low 
charged particle multiplicity and satisfies the tau cri- 
teria. A small additional contribution to the back- 
ground comes from W and Z events where the tau 
jet comes from a tau or mismeasured electron from 



the boson decay. In this case the tau jet typically 
has one or three associated charged tracks. 

A combination of events satisfying the various 
jet energy triggers in the experiment models the 
hadronic background well. The background normal- 
ization is computed function of the Et and 
charged multiplicity of the tau. The normaliza- 
tion equalizes the number of events of any charged 
multiplicity except 1 or 3, in three ranges of Et- 
The Monte Carlo simulation shows that real taus 
contribute less than a few percent to these bins in 
multiplicity. An estimated total of 17.4ib2.5 events 
come from processes where the tau jet came from an 
hadronic jet; the error is statistical only. 

The estimate of the non-hadronic-jet contribution 
to the background comes from Monte Carlo simula- 
tion of the various processes. Of these, only the con- 
tribution from Z ^ T'^T~ remains non- negligible af- 
ter all cuts. Using a total of 30, 000 events generated 
with the ISAJET program and then passed through 
the CDF detector simulation and analysis, we expect 
l.libO.4 (stat) events with 1 or 3 associated charged 
particles. The production cross section comes from 
the measured Z cross section, assuming lepton uni- 
versality: a{pp ZX;Z^e+e-X) =0.209 ± 
0.013(stat) ± 0.017(syst) nb g. This background 
is small for several reasons: the process has a small 
cross-section, the two outgoing taus are azimuthally 
back-to-back, and the is typically not large since 
the tau decay neutrinos are back-to-back. 

The predicted background from W-l-jets events, 
with W Tu, comes from 40,000 HERWIG events 
which were passed through the CDF detector sim- 
ulation, including the relevant trigger efficiencies. 
Again assuming lepton universality the produc- 
tion cross section a{pp^WX;W^ev) = 2.19± 
0.04(stat)±0.21 (syst) nb |jll| is used for normaliza- 
tion. Most of these events are rejected by the cuts 
on S{^t), and No event passed the se- 

lection criteria. 

The other processes involving W and Z bosons 
result in background taken into account by the 
hadronic jet sample, and contribute negligibly to the 
non-hadronic-jet component. 

There is a small acceptance for events from stan- 
dard model top quark pair production; for a top 
quark with mass of 176 GeV/cP one expects 0.2ib0.1 
(stat) events. This acceptance affects the number of 
expected events in the signal, and does not enter the 
background estimate. 
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Figure 4: Charged particle multiplicity distribution 
in the 10° cone after all the cuts, for the data sam- 
ple (points), for the background normalized to the 
data (cross-hatched histogram), and for the expected 
signal from the normalized Monte Carlo simulations 
(open histogram) added to the background. 

VII. SYSTEMATIC UNCERTAINTIES 

Systematic effects which can lead to uncertainty 
in the final result can be classified into those which 
affect the background estimate and those which af- 
fect the number of expected events. Many of the 
systematic uncertainties affecting the number of ex- 
pected events depend on the top quark mass. Table 
II lists the different estimated systematic uncertain- 
ties. For the cases where there exists a top quark 
mass dependence, the extreme values appear in the 
table. 

Various effects can bias the background estima- 
tion, such as the binning of the Et distribution and 
the normalization method. Dividing the Et distri- 
bution into smaller bins and following the same nor- 
malization method leads to a negligibly small differ- 
ence in the expected number of hadronic background 
events. The normalization method is based on jet 
trigger events, but one can check for a trigger bias 
by removing the jet which was responsible for the 
trigger, and no significant effect appears. The to- 
tal hadronic background is conservatively estimated 
to be 17.4ib2.5(stat)±0.6(sys) events, based on these 
cross checks. 

We have compared the number of expected events 
for different masses of the top quark and the charged 



Top quark cross section 30-10% 

ISAJET gluon radiation 16-1.3% 
Integrated luminosity 3.7% 
Trigger efficiency 5.5% 
MC statistics 10-4% 

Energy scale effect 32-7.5% 
Background estimation 14% 

Table II: Sources and magnitudes of the systematic 
uncertainties in the analysis. The values are the rel- 
ative uncertainties in the number of expected events, 
and represent the extremes for the top mass range 
100-174 GeV/c2. 



Higgs with and without initial-state gluon radiation 
in ISAJET. Half the difference between these num- 
bers was taken as the systematic uncertainty. The 
mean value between the number of expected events 
with and without initial-state gluon radiation was 
taken as the number of expected events. The isola- 
tion cut is the criterion most affected by initial-state 
gluon radiation; the number of jets is smaller with 
no initial-state gluon radiation. The probability to 
have associated charged particles of the tau candi- 
date mixed with a particle of another jet becomes 
smaller and the efficiency of the cut increases. This 
effect also depends on the top quark mass, since for 
a heavy top quark there is less energy to produce 
jets as a result of the initial-state gluon radiation. 

The systematic uncertainty on the trigger ef- 
ficiency was estimated by varying each point of 
the measured trigger efficiency by its uncertainty. 
The relative uncertainty on the number of expected 
events due to the systematic uncertainty in the trig- 
ger efficiency is conservatively estimated to be 5.5%. 
In this calculation, mtop = 120 GeV/c^ and mniggs = 
100 GeV/c2. 

The absolute energy scale uncertainty varies from 
±10% at 8 GeV to ±3% at 100 GeV. In the Monte 
Carlo simulations, we shifted the jet energy scale by 
these values, and repeated the analysis, reconstruct- 
ing the Et of each jet and other relevant event pa- 
rameters. We used the mean relative difference of 
the change in the number of expected events when 
the energy scale is shifted as the systematic uncer- 
tainty on the energy scale. 
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VIII. RESULTS 



After selection, there remain 74 events from the 
data sample, of which a total of 19 events have a tau 
candidate with either 1 or 3 associated charged par- 
ticles. Fig. ^ shows the multiplicity distribution for 
the data sample and the hadronic background nor- 
malized to the data. For comparison the plot shows 
the distribution from the Monte Carlo signal simu- 
lation normalized to the total integrated luminosity 
and added to the hadronic background estimation. 
The estimated total number of background events 
is 18. 5 lb 2. 6, where the error comes from adding in 
quadrature the systematic and statistical uncertain- 
ties. 

The mass limits must take into account the un- 
certainties, both statistical and systematic, on the 
number of expected background and signal events. 
For a given mass point, a simple Monte Carlo gen- 
erates a large ensemble of trials with the numbers of 
expected signal and background events varying in a 
Gaussian fashion about the mean. In each trial it 
generated a number of observed events from a Pois- 
son distribution with a mean equal to the number 
of signal plus background. The standard deviations 
of the Gaussians are the combined statistical and 
systematic uncertainties. The mass point can be ex- 
cluded with 95% confidence if in 95% or more of 
the trials the total number of events exceeds the 19 
events actually observed. 

Fig. H shows the resulting limit for large values of 
the parameter tan /3, for which the branching ratios 
of t — > Hb and H ^ tv approach unity. Fig. |^ shows 
the limit for tan (3 = 50, 100, and 500 in the plane of 
the top quark mass versus the charged Higgs mass. 
As the charged Higgs mass decreases, the missing 
transverse energy decreases on average, reducing the 
efficiency of the selection. Also, as the top mass 
increases, its production cross section decreases, and 
the number of expected events decreases. As the 
parameter tan (3 increases, the branching ratios of 
the top to charged Higgs and charged Higgs to tau 
both increase, allowing a better limit. 

The event selection used is well optimized for large 
masses of the top quark and the charged Higgs. The 
present statistics exclude the region for large tan /3, 
extending from the limit of the previous analyses 
just up to the region where the top mass has been 
measured. 
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Figure 5: Regions of the (w-topj^-Higgs) plane ex- 
cluded at 95% CL for a 100% branching ratio of 
— > TV. The plot also shows the limit from the 
previous analyses [Q, ^. 
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